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Abstract 

We argue that the recently introduced "statefinder parameters" (Sahni et al., JETP Lett. 77, 201 
(2003)), that include the third derivative of the cosmic scale factor, are useful tools to characterize 
interacting quintessence models. We specify the statefinder parameters for two classes of models 
that solve, or at least alleviate, the coincidence problem. 
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I. INTRODUCTION 



In the last years the conviction that our Universe is undergoing a phase of accelerated 
expansion has gained further ground among cosmologists albeit the nature of the cosmic 
substratum (usually called dark energy) behind this acceleration remains as elusive as ever 
P]. While several candidates have been proposed ^ there is no agreement on which of them 
should be considered as the favored one. By all accounts, much more observational input is 
needed before this point might be settled. 

The acceleration is evaluated by the deceleration parameter q = — a/(aiJ^), where a{t) 
stands for the scale factor of the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 
and H = a /a for the Hubble parameter. As mentioned above, the present state of the 
Universe seems to be characterized by a negative q, but it is hard to determine its value 
observationally. Therefore, since many models predict acceleration some further information 
should be welcome. 

Among others, there are cosmological models whose evolution is dominated by in- 
teracting components -say, dark matter and dark energy. Models in which the main 
energy components do not evolve separately but interact with each other bear a special 
interest since they may alleviate or even solve the "coincidence problem" that afflicts many 
approaches to late acceleration P]. This problem can be summarily stated as "why now?", 
that is to say: "Why the energy densities of the two main components happen to be of 
the same order today?" In this paper we focus on a Universe filled with two components, 
namely, non-relativistic matter (subscript m) with negligible pressure, i.e., Pm <^ Pm and 
dark energy (subscript x) -with equation of state = wp^ where w < -, such that the 
latter decays into the former component according to 



Pm + 3Hpm = Q , 

p^ + 3H{i + w)p^ = -g, (1) 

where Q > measures the strength of the interaction. For later convenience we will write 
it as Q = —3IIH where the new quantity 11 has the dimension of a pressure. 
The Einstein field equations for spatially flat FLRW cosmologies are 
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where p = + is the total energy density, and we have set c = 1. The quantity H 
is related to the deceleration parameter g by g = — 1 — {H / H'^) = (1 + 3wQx)/'^, where 
= SnCpx/iSH"^). It is obvious, that the deceleration parameter does not depend on 
whether or not both components are interacting. However, differentiating H again, we 
obtain 



H 9 / pA 9 



Px 



n 



w 



w{l + w} w— 

p p 3H p 



(4) 



At variance with H and H, the second derivative H does depend on the interaction between 
the components. Consequently, to discriminate between models with different interactions 
or between interacting and non-interacting models it is desirable to characterize the cosmo- 
logical dy,>am,cs additionally by parameters that depend on H. Recently. Sahni et al. Q 
and Alam et al. p have introduced a pair of new cosmological parameters (the so-called 
"statefinder parameters") that seem to be promising candidates for this purpose. These are: 
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In the present context of interacting fluids they take the form 
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where k, = Pm/Px, and 
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For non-interacting models i.e., for 11 = 0, these parameters reduce to the expressions studied 



3 



m and 0. The target of this short communication is to illustrate how the statefinder 
parameters may be of help when exploring cosmological models whose evolution is dominated 
by interacting components. While we focus ourselves on interacting cosmologies, we mention 
that the third derivative of the scale factor is generally necessary to characterize any variation 
in the overall equation of state of the cosmic medium This becomes obvious from the 
general relation 
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where P is the total pressure of the cosmic medium which in our case reduces to P ^ px- 
Since 
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it is evident, that an interaction term in P ^ px = wpx + wpx according to (P) will addi- 
tionally change the time dependence of the overall equation of state parameter P/ p. 

Interacting models allow a dynamical approach to the coincidence problem. The central 
quantity here is the density ratio k introduced beneath Eq. This parameter should 

be a constant of the order of unity at late times for the coincidence problem to be strictly 
solved. At least it should, however, vary slowly over a time of the order of H~^. The ratio 
K is governed by the evolution equation (cf. Eqs. (^) 
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Below we study the statefinder parameters for different solutions of this equation, corre- 
sponding to two broad classes of matter-quintessence interacting models. 

II. SCALING SOLUTIONS 



In a recent paper ^ the authors showed that scaling solutions, i.e., solutions of the form 
Pm/ px oc a~'^, where ^ denotes a constant parameter in the range [0,3], can be obtained 
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when the dark energy component decays into the pressureless matter fluid -Eqs. (jT]). These 
solutions are interesting because they alleviate the coincidence problem . Indeed, a model 
with ^ = 3 amounts to the ACDM model with w = —1 and 11 = 0. For the opposite extreme 
value ^ = the Universe dynamics admits a stable, stationary solution k = constant, 
thereby no coincidence problem arises Q]. Hence, the deviation of the parameter ^ from 
^ = quantifies the severity of the problem whereby any solution that differs from ^ = —3w 
represents a testable, non-standard cosmological model and any solution with < 3 renders 
the coincidence problem less acute. In that scheme, with w = constant, it can be shown 
that the interactions which produce scaling solutions are given by 

where denotes the present energy density ratio and z = (ao/a) — 1 is the redshift. Inserting 
this expression in Eqs. © and ((Tj) we get for the statefinder parameters 
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Figure 1 depicts the function r(s) for different values of ^. The lower the value of ^, the 
lower the value of the corresponding curve in the s-r plane and the less acute results the 
coincidence problem. Note that qualitatively these curves remain similar to those of non- 
interacting models (see Fig.l of Ref.|5j]). For the sake of comparsion, we note that the ACDM 
model (n = 0, If = —1) corresponds to the point (not shown in the figure) s = 0, r = 1. 
From an observational point of view and for discriminating between different models, the 
current values tq and Sq of the statefinder parameters are of particular relevance. For the 
scaling models we have 

9 w , ^ ( i\ /, .\ 

^^0 = 1+0^"^ ^0, and sq = 1 + w - \w + -\— . (14) 
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FIG. 1: Selected curves r(,s) in the redshift interval [0,6] (from left to right) with w = —0.95 and 
kq = 3/7 for three different values of the parameter ^, viz (a) 2.5; (b) 1.5; (c) 0.5. 



Realizing that 
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and introducing 



qoA = qo{w^ -1) = - 



12 — kq 3 kq 



21 + kq 21 + kq 



1 + ?0A , 



(16) 



we may classify the different models through their dependence so(?o) which reads 



go - goA + ( 2 - 1 ) (1 + ^oa) 



(17) 



The first part in the bracket on the right hand side describes the deviation from w = —1, 
the second part accounts for the deviations from the ACDM scaling ^ = 3. For models 
with w — —1, e.g., which all have the same deceleration parameter go = Qqa, we have 
So = § (I — l) (1 + ?oa)- Of course, ^ = 3 corresponds to the ACDM model with sq = 0. 
Assuming kq — 3/7, a scaling ^ = 1 results in sq (^ = 1) = —0.2, while the stationary solution 
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^ = has So = 0) = —0.3. Similar considerations hold for other values of w. Thus, the 
parameter sq is able to discriminate between different scaling models, characterized by the 
same deceleration parameter. 



A. Luminosity distance 

It is interesting to see how the statefinder parameters enter the expression for the lumi- 
nosity distance. Up to second order in the redshift the Hubble rate is 



H{z)=Ho+ [ — ] z + 
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By virtue of 
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this can be written as 
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The luminosity distance 
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dz 



becomes (cf. y\) 
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Since the interaction affects ro but neither go nor Hq, it is obvious that the luminosity dis- 
tances of different interacting as well as of interacting and non-interacting models manifests 
itself only in third order in the redshift. For the scaling model leading to ^ = 1, and w = —1 



7 



we have go = — (2 — /to)/[2(H-Ko)] and tq = 1 — [3ko/(1 + ft;o)^], and we recover the previously 
obtained expression (see Eq. (38) of Ref. jsl) 
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For the ACDM model the factor (6 + kq) occurring in last expression is replaced by (10 + kq). 
On the other hand, for a given w all the scaling models, including ACDM, are degenerate 
with respect to the deceleration parameter. For all these models the present value of q is 
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This expression also holds true for the asymptotically stable model of the next section. This 
demonstrates explicitly that discrimination between interacting and non-interacting models 
or between different interacting models requires the knowledge of the luminosity distance 
up to the third order in the redshift. In other words, the circumstance that the interaction 
is felt only by parameters containing the third derivative of the scale factor corresponds to 
the fact that the luminosity distance of interacting models is affected only in the third order 
of the redshift. 



III. ASYMPTOTICALLY STABLE SOLUTIONS 

For the special case that the interaction term is assumed to obey IT = — c^p with (? = 
constant < 1, the evolution equation ()10|) has two stationary solutions for w = constant, 
namely, and = l/nf (with > 1), given by 



2c2 



w 



=F \/w{w + 4c2) 



(25) 



-see Ref. [ll| for details. It can be shown that whereas the solution is stable the solution 
nt is unstable. There exists a solution 
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with y — (aeq/a)^ where Ueq is the scale factor at which the energy density of dark energy 
equals the energy density of dark matter and A — — 3w(l — k,j)/{1 + k,'^), according to which 
K evolves from a matter dominated phase (k+ > 1) for a ^ a^q to a dark energy dominated 
phase {kJ < 1) for a ^ a^q as the Universe expands. In this case we have 



(27) 



with 7/0 = «^(«o - i^s )/{'^t - f^o)- 

Accordingly, we readily obtain that the statefinder parameters of this model reduce to 
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and 
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where (? = —wn~/ (1 + is valid. 

Figure 2 shows some graphs of the function r(s) for different choices of k~ . At variance 
with graphs of the scaling model, the location of the curves in the plane (s, r) is unrelated 
to the alleviation or solution of the coincidence problem since k — kJ — constant at late 
times for all the cases of this model. 

For the present value of the parameter s this model yields 



So = 1 + + (1 + 'Co) , where 1 + Kq — {1 + k^) 
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(30) 



A parallel study to the case of the scaling solutions leads to 



So = 



(31) 



Again, the first two terms in the bracket account for the difference to models with w = —1 
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FIG. 2: Selected curves r(s) in the redshift interval [0,6] (from left to right) with w = —0.95 and 
kq = 3/7, for different values of kJ, viz (a) 0.3; (b) 0.35; (c) 0.4. 

for which go = Qoiw = —1) = Qqa- The term describes the impact of the interaction on 
the parameter sq. 

IV. CONCLUDING REMARKS 

n ■-' 

The statefinder parameters introduced in |3| and |6| are expected to be useful tools in 
testing interacting cosmologies that solve or at least alleviate the coincidence problem which 
besets many approaches to late acceleration. It is manifest that while the deceleration pa- 
rameter does not feel the interaction between the dark energy and dark matter the statefinder 
pair (r, s) does. 

We hope that in some not distant future we will have at our disposal observational 
techniques capable of determining these parameters. These are bound to shed light on the 
nature of dark energy and dark matter. 
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